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SUMMARY 


The  flow  of  a  compressible  fluid  through  a  channel  having  locally 
supersonic  regions  is  studied  by  using  the  Tricomi  equation  in  the  hodo- 
graph  variables  as  an  approximation  in  the  sonic  region  to  the  equation 
of  flov  of  an  irrotational,  inviscid  gas.  It  is  shown  that  this  is 
equivalent  to  studying  the  flow  of  a  gas  having  a  pressure-density  rela¬ 
tion  matching  the  isentropic  relation  to  the  third  derivative  at  the 
sonic  point.  A  one-parameter  family  of  solutions  of  the  Tricomi  equa¬ 
tion  is  used  which  provides  symmetrical  accelerated-decelerated  flows. 
The  variation  of  this  parameter  alters  the  Mach  number  at  the  center  of 
the  throat,  the  velocity  distribution  and  gradient  along  the  center 
streamline,  as  well  as  the  shape  of  the  channel,  that  is,  the  curvature 
of  the  bounding  streamline. 

As  specific  examples,  flows  are  computed  having  Mach  number  equal 
to  unity  and  to  0.86  at  the  center  of  the  throat  section.  Constant- 
velocity  lines  are  plotted  and  it  is  found  that  the  velocity  gradient 
becomes  zero  at  three  places  along  each  streamline  outside  of  a  limiting 
streamline  for  values  of  the  parameter  greater  than  zero  (M<1  at 
center  of  throat  section).  For  the  parameter  equal  to  zero  (M  =  1  at  > 
center  of  throat  section) ,  the  velocity  gradient  along  the  streamlines 
and  the  curvature  is  discontinuous  at  all  points  of  the  two  character¬ 
istics  which  meet  the  center  streamline. 

Other  solutions  to  the  Tricomi  equation  are  discussed  which  may  be 
used  to  formulate  channel  flows.  The  exact  nature  of  these  flows  has 
not  yet  been  investigated. 


INTRODUCTION  ■ 

In  the  investigations  of  the  effect  of  free-stream  Mach  number  on 
the  lift  and  drag  of  an  airfoil  in  a  uniform  stream,  it  was  found  that 
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the  lift  did  not  begin  to  decrease  or  the  drag  increase  until  the  free- 
stream  velocity  exceeded  the  critical  free-stream  velocity,  that  is, 
the  velocity  for  which  the  local  velocity  of  sound  is  just  attained  at 
some  point  of  the  body.  Thus  it  appeared  that,  for  free-stream  velocities 
in  between  the  critical  free-stream  velocity  and  the  velocity  for  which 
the  lift  characteristics  changed  rapidly  with  Mach  number,  there  might 
exist  about  the  body  a  locally  supersonic  region  which  contained  no  shocks, 
although  such  shock- free  flows  have  never  been  observed  experimentally. 

A  similar  phenomenon  occurs  in  the  study  of  flows  through  a  gradual 
constriction.  When  the  exit  pressure  of  a  nozzle  is  decreased,  the 
mean  velocity  at  the  narrowest  cross  section  and  also  the  total  mass 
flow  are  increased.  However,  when  the  mean  velocity  at  the  constriction 
reaches  the  local  velocity  of  sound,  the  mass  flow  reaches  a  maximum; 
and  the  flow  then  changes  from  an  accelerated-decelerated  flow  to  one 
in  which  the  fluid  is  always  accelerated  in  passing  through  the  nozzle. 
Further  decreases  in  the  exit  pressure  fail  to  increase  the  mass  flow 
but  only  introduce  shocks  into  the  supersonic  side  of  the  nozzle  to 
satisfy  the  conditions  of  conservation  of  mass,  momentum,  and  energy. 

One  question  which  occurs  regarding  such  a  process  is:  Do  there  occur 
isolated  regions  of  shock-free  supersonic  flow  in  the  neighborhood  of 
the  walls  before  the  mass  flow  reaches  the  maximum  and  the  flow  becomes 
purely  accelerated? 

This  work  was  done  at  Brown  University  under  the  sponsorship  and 
with  the  financial  assistance  of  the  National  Advisory  Committee  for 
Aeronautics. 


SYMBOLS 

x,y  Cartesian  coordinates  (z  =  x  +  iy) 

u,v  velocity  components  in  x-  and  y-direction,  respectively 

7  ratio  of  specific  heats 


a^_  Mach  angle 

t  =  q^/23  (t  =  1  corresponds  to  limiting  speed) 

q  velocity  vector 

q  magnitude  of  velocity  vector 


* 
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© 

yi 


0  inclination  of  velocity  vector  with  x-axis 

T-l  =  2t(  1  -  T)"P 


t2 


1  -  (23  +  1)t 
2t( 1  -  t)P+1 


c 


p 

V 

p 

Ps 

M 


local  speed  of  sound  normalized  to  speed  of  sound  in 
still  air 

density  of  fluid  normalized  to  stagnation  density 
specific  volume 

pressure  of  fluid  corresponding  to  density  p 
critical  pressure  (pressure  at  s  =  0) 

Mach  number  (q/c) 

2,1120s* 

1  +  1.2304s* 


velocity  potential 
stream  function 

(1.45320  -  0.84664s)t 

conjugate  potential  to  (by  equations  (6)) 

symmetrical  stream  function  given  "by  equation  (19) 
stream  function  given  by  equation  (15) 


9 

♦  " 

\|rf  =  = 

9* 

f 

♦l 
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EARLY  INVESTIGATIONS  OF  LOCALLY  SUPERSONIC  REGIONS 


The  first  attempt  to  construct  a  flow  through  a  channel  having 
locally  supersonic  regions  was  made  hy  Taylor  in  1931  (reference  1). 

He  expanded  the  velocity  potential  in  a  double  power  series  about  some 
central  point  of  a  two-dimensional  nozzle  whose  walls  consisted  of  two 
arcs  of  circles.  Taylor  was  able  to  obtain  several  flows  which  were 
symmetric  in  the- velocity  about  the  line  of  centers  of  the  two  circles. 
The  values  of  the  maximum  velocity  on  the  straight  streamline  for  these 
flows  were  bounded  by  a  certain  velocity  for  which  the  local  velocity 
of  sound  on  the  wall  wos  just  barely  exceeded.  From  this  result, 
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Taylor  concluded  that  the  potential  flow  through  a  channel  with  locally 

*  supersonic  regions  was  impossible.  Gortler  (reference  2) ,  however, 
felt  that  Taylor's  conclusion  was  unfounded,  for  if  more  terms  of  the 
power  series  were  taken  and  the  boundary  conditions-  on  the  wall  satisfied 
to  more  terms  in  the  series  expansion,  then  the  difficulty  which  Taylor 
encountered  in  determining  the  coefficients  might  not  exist.  Gortler, 
using  much  of-  Taylor's  computations,  satisfied  one  less  relation  of  the 
boundary  conditions  than  Taylor  and  thereby  had  two  free  parameters  to 
vary  instead  of  one.  One  of  these  parameters  represented  the  velocity 

at  the  center  axis  of  the  throat  Of  the  nozzle,  and  the  other,  the 
velocity  gradient  along  the  center  axis  at  this  point.  By  varying  the 
velocity  gradient  Gortler  obtained  a  sequence  of  channels  which  exhibited 
locally  supersonic  regions  near  the  channel  walls  at  the  constriction. 

As  the  velocity  gradient  was  increased,  the  locally  supersonic  region 
gradually  extended  toward  the  axis.  When  the  local  velocity  of  sound 
was  attained  at  the  axis,  the  flow  became  identical  to  Taylor's  unsym- 
metrical  flow  through  a  Laval  nozzle.  Gortler  showed  that  the  locally 
supersonic  region  cannot  extend  to  the  center  axis  if  the  velocity 
gradient  is  to  be  continuous  along  the  center  streamline.  When  the 
s  sonic  velocity  is  just  attained  on  the  axis,  but  not  exceeded,  then  the 

condition  of  continuity  of  the  gradient  will  require  a  zero  gradient  at 
^  the  sonic  point.  Under  this  condition  the  following  theorem,  proved  by 

Gortler,  shows  that  the  sonic  line  must  be  a  straight  line  and  the 
curvature  of  the  wall  at  the  throat  must  be  zero. 

Theorem;  Let  an  inviscid  flow  of  a  compressible  gas 
be  given  which  is  symmetrical  about  the  y  =  0  axis  and  is 
analytic  in  the  neighborhood  of  a  point  x  =  x0  and  y  =  0. 

Further,  suppose  that  the  local  velocity  of  sound  is  attained 
at  this  point  of  the  center  axis  and  that  the  velocity  gra¬ 
dient  vanishes;  that  is,  u(xD,o)  =  1  and  |^(x0,o)  =  0. 

Then  it  follows  necessarily  for  all  y1 s  in  the  neighbor- 
hood  that  u(xQ,y)  =  1,  |j(xQ,y)  =  0,  and  v(x0,y)  s  0. 

Flows  having  locally  supersonic  regions  in  which  the  sonic  velocity  is 
just  reached  on  the  center  axis,  therefore,  will  not  be  analytic  in  the 
neighborhood  of  the  center-axis  point  of  the  throat.  Such  a  flow  has 
been  computed  by  Tomotika  and  Tamada  and  a  similar  flow  having  symmetry 
about  both  horizontal  and  vertical  axes  is  given  in  this  work.  The 
velocity  gradient  in  both  flows  has  a  jump  at  the  center  of  the  nozzle. 

The  work  of  Tomotika  and  Tamada  (reference  3)r  using  an  approximate 

*  gas  law,  and  of  Frankl  (reference  4),  Gortler  (reference  2),.  Oswatitsch 
and  Roths te in  (reference  5),  Fox  and  Southwell  (reference  6),  and 
others,  using  power  series,  already  indicates  strongly  that  solutions 
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which  have  locally  supersonic  regions  for  the  exact  flow  equations  do 
exist.  Recognition  must  be  given,  of  course,  to  the  lack  of  knowledge 
about  the  convergence  of  the  power  series,  and  to  the  fact  that  the 
Tricomi  equation  which  the  Japanese  a.uthors  use  approximates  the  actual 
gas  only  in  a  small  neighborhood  of  the  sonic  velocity.  Yet,  the  fact 
that  both  methods  give  flows  having  locally  supersonic  regions  clearly 
points  toward  a  similar  result  from  the  exact  equation  for  an  isentropic 
compressible,  inviscid  fluid.  But  the  existence  of  these  solutions  by 
no  means  settles  the  matter.  It  is  to  be  noted  that  "existence"  is 
used  in  the  mathematical  sense.  Actually,  all  experimental  results 
(Liepmann  (reference  7)  and  Frankl  (reference  4))  show  that  irrotational 
flow  breaks  down  soon  after  the  appearance  of  the  lower  critical  Mach 
number  (Tsien  and  Kuo  (reference  8)).  None  of  the  experiments  shows 
locally  supersonic  regions  of  any  appreciable  extent  without  shocks. 


GENERAL  PROPERTIES  OF  A  LOCALLY  SUPERSONIC  REGION 


Frankl  /reference  9)  and  Nikolsky  and  Taganoff  (reference  10)  have 
made  some  general  investigations  of  the  properties  of  the  locally  super¬ 
sonic  regions  without  shocks.  Frankl  (reference  9)  showed  that  the 
symmetric  flow  over  an  obstacle  which  is  symmetric  about  two  perpen¬ 
dicular  axes  is  uniquely  determined  by  the  boundary  conditions  on  the 
portions  of  the  arcs  DFD!  and  EGE*  in  figure  1  which  lie  outside  of  the 
region  defined  by  the  two  pairs  of  characteristics  BD  and  BE  and  B* D* 
and  BTET  and  the  conditions  at  infinity.  The  segments  of  the  body  DE 
and  D'E1  between  the  characteristics  then  depend  upon  the  rest  of  the 
body  and  the  free-stream  Mach  number.  Thus  it  appears  that  symmetric 
flow  with  supersonic  regions  does  not  exist  for  an  arbitrary  body, 
although  for  a  given  free-stream  Mach  number  there  might  well  be  a 
large  class  of  profiles  for  which  such  flows  do  exist. 

The  results  of  Frankl  can  be  shown  more  rigorously  and  easily  by 
making  use  of  the  following  uniqueness  theorem  proved  by  Tricomi 
(reference  ll): 

In  a  domain  bounded  by  a  curve  in  the  subsonic  region 
with  its  termini  on  the  sonic  line  in  the  hodograph  plane 
and  two  intersecting  characteristics  drawn  from  these  ter¬ 
mini,  the  solution  of  the  Tricomi  equation  is  uniquely  deter¬ 
mined  by  the  values  prescribed  on  the  subsonic  arc  and  on 
one  of  the  characteristic  segments  bounding  the  domain. 

The  character  of  the  flow  through  a  nozzle  having  locally  supersonic 
regions  is  shown  in  the  hodograph  plane  in  figure  2.  If  the  value  of 
the  stream  function  is  given  on  the  arc  BCD  and  the  segment  of  the 
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axis  OAB  and  the  characteristic  ED,  then  the  solution  is  given  uniquely 
inside  the  region  BCDEOB.  In  the  physical  plane,  this  means  that  the 
solution  is  given  uniquely  in  the  region  A'B’C'D'E'O' A'  (fig.  3)  since 
the  arc  OA  in  figure  2  maps  into  the  line  9  =  0  (O' A'  in  fig.  3)  the 

line  AB  into  the  center  streamline,  and  the  characteristics  DE’and  EO 
into  the %Mach  lines  D'E*  and  E'O’.  The  same  thing  can  he  shown  for 
A  0  G'H'I'J'A'  since  symmetry  is  assumed.  Since  the  values  of  ©,  q, 
and  \|r  are  known  on  the  two  Mach  lines  E'O'  and  O'G' ,  then  the  solution 
is  determined  inside  the  quadrilateral  E’L'G'O’.  Thus  the  flow  field 
is  entirely  determined  by  the  boundary  conditions  on  the  arcs  C*F' 
and  M»I',  and  by  the  velocity  distribution  on  the  center  streamline  JB' 
and  on  C'B'  and  I'J'. 

The  fact  that,  for  an  inviscid  flow  which  has  locally  supersonic 
regions,  the  boundary  cannot  be  chosen  in  an  arbitrary  manner  actually 
has  little  significance  in  the  flow  of  real  gases.  In  such  flows,  there 
exists  a  boundary  layer  near  the  walls  and  a  pertinent  question  regarding 
locally  supersonic  regions  may  be  stated  as  follows;  Does  the  boundary 
layer  adapt  itself  so  that  locally  supersonic  regions  can  occur  without 
shocks,  or  are  shocks  always  present  in  locally  supersonic  regions? 

Nikolsky  and  Taganoff  have  discovered  some  of  the  properties  which 
a  locally  supersonic  region  in  an  inviscid  flow  free  of  shock  waves 
must  have.  Briefly  these  properties  are  as  follows: 

(1)  The  slope  of  the  contour  at  any  given  point  in  the  locally 
supersonic  region  is  the  arithmetic  mean  of  the  slopes  of  the  velocity 
vector  at  the  points  on  the  transition  line  which  are  the  termini  of 
the  characteristics  ’originating  from  the  given  point  of  the  contour. 

(2)  If  one  moves  along  the  transition  line  in  such  a  way  that  the 
region  of  subsonic  velocities  lies  on  the  left,  then  the  velocity  vector 
will  always  rotate  clockwise. 

(3)  It;  in  the  supersonic  zone,  a  segment  of*  a  characteristic  of 
one  family  is  given  such  that  the  characteristics  of  the  other  family 
originating  from  points  of  this  segment  extend  to  the  transition  line, 
then  the  inclination  of  the  velocity  vector  and  its  magnitude  are 
monotonic  functions  along  the  given  segment  of  the  characteristic. 

(4)  The  magnitude  of  the  velocity  along  a  rectilinear  section  of 
the  profile  in  a  supersonic  zone  must  decrease  in  the  direction  of  flow, 

(5)  Prandtl-Meyer  flow  between  two  characteristics  of  one  family, 
with  the  characteristics  of  the  other  family  rectilinear,  cannot  be 
completely  realized  up  to  the  transition  line  in  the  finite  region. 
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(6)  Characteristics  of  the  first  family  which  originate  from  points 
of  a  rectilinear  section  of  the  profile  cannot  extend  to  the  transition 
line,  but  must  instead  terminate  in  a  compression  shock.  * 

Nikolsky  and  Taganoff  also  discuss  the  possibility  of  the  Mach  lines 
forming  an  envelope  in  the  supersonic  zone.  They  show  that  if  the  solu¬ 
tion  to  the  flow  over  a  body  is  found  which  satisfies  the  boundary  con¬ 
ditions  then  the  Jacobian  of  the  transformation  from  the  hodograph  plane 
to  the  physical  plane  cannot  vanish  on  the  boundary,  since  the  curvature 
of  the  streamline  must  be  zero  when  J  =  0.  By  using  the  property  of 
statement  (3)  above  they  show  that  if  the  Jacobian  does  not  vanish  on 
the  boundary  then  it  cannot  vanish  in  the  locally  supersonic  region  in 
the  flow  field  outside  the  body.  They  give  as  the  criterion  of  break¬ 
down  the  relation 


—  =  q  tan  a*, 
&Q  ± 


where  q  is  the  velocity  magnitude,  G  the  angle  of  inclination  of  the 
contour,  and  cx]_  the  Mach  angle  associated  with  the  local  Mach  number. 

What  Nikolsky  and  Taganoff  failed  to  notice,  however,  is  that  the  above 
equation  is  the  condition  that  the  streamline  representing  the  body  be 
tangent  to  a  characteristic  in  the  hodograph  plane.  If  the  derivatives 
with  respect  to  6  and  q  of  the  stream  function  are  finite,  then  the 
streamline  will  have  infinite  curvature  at  the  point  where  the  relation 
above  is  satisfied.  What  Nikolsky  and  Taganoff  have  actually  shown  may 
be  stated  more  precisely  in  the  following  form: 

Theorem:  If,  along  the  portion  of  a  streamline  lying  in 

a  local  supersonic  region  of  flow,  there  exists  no  point  at 
which  the  curvature  becomes  infinite,  then  there  cannot  exist 
any  limit  line  inside  the  supersonic  zone  lying  between  the 
given  streamline  and  that  streamline  which  is  just  tangent 
to  the  transition  line. 

However,  if  one  considers  the  continuation  of  the  flow  pattern  then  one 
may  find  a  streamline  for  which  the  curvature  is  infinite  and  hence  the 
Jacobian  J  =  0.  Nikolsky  and  Taganoff  were  unable  to  show  that  J  =  0 
on  the  transition  line  was  an  impossibility  under  the  conditions  of  the 
theorem  above.  Friedrichs  (reference  12),  however,  showed  that,  under 
the  hypotheses  of  the  theorem  above,  infinite  acceleration  and  infinite 
curvature  of  the  streamlines  could  not  occur  on  the  transition  line 
between  the  two  streamlines  mentioned  in  the  theorem.  In  investigating 
the  proximity  to  breakdown  of  a  potential  flow  over  an  obstacle,  it 
seems  advisable  to  study .the  continuation  of  the  flow  into  the  region 
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inside  the  body  and  investigate  how  near  to  the  body  streamline  lie  the 
cusps  of  a  limit  line.  If  limit  lines  lie  near  the  body>  then  it  is 
probable  that  a  slight  change  in  the  boundary  or  an  increase  in  the 
free  stream  would  result  in  an  inability  to  satisfy  the  boundary 
conditions. 


APPROXIMATE  HODOGRAPH  METHOD  FOR  TRANSONIC  FLOW 
INVOLVING  SOLUTIONS  TO  THE  TRICOMI  EQUATION 


In  this  discussion  an  approximation  to  the  compressible-flow  hodo- 
graph  equations  will  be  used  in  which  the  Tricomi  equation  replaces  the 
usual  second-order  linear  equation  for  the  stream  function.  In  the 
notation  of  reference  13,  the  equations 


<Pe  =  Tpt. 


cpT  =  -T2t0 


where  Tl  -  2t(1  -  t)"3  and  T2  =  — ..(2P  +  lK  =  represent  an 

2t( 1  -  t)^+^  p 

exact  linearization  of  the  equations  of  irrotational  motion  of  an 
inviscid,  compressible  fluid.  If  a  new  variable 


(20+1) 


(1  -  t)P 


is  introduced  into  equation  (  l)  ^  the  equations 


*0  ■  "V 


<P.»  -  0)*(s*He 
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are  obtained, 
equation  for 


where  co*(s*)  = 
\|r  results  when 


W 


Then  the  following  differential 


cp  is  eliminated  from  equation  (2): 


^s*s*  +  =  0 


(3) 


The  general  forms  for  two  simultaneous  differential  equations  in 
two  independent  variables  which  reduce  to  the  Tricomi  equation  have  been 
developed  by  Loewner  (reference  l4) .  One  such  form  will  be  considered 
here.  Let  \|r  and  cp  be  defined  by  the  relations 


Vs  =  4  ^ 


cps  =  sco^e 


(*) 


and  let  a  transformation  of  the  dependent  variables  be  given  by 


\| t'  =  a 


'P’s  =  -~VS 


.  1  dah 

cp’e  =  -  —  9©  +  —  * 
u  °  ds 


J 


(5) 


On  eliminating  cp*  from  equation  (5),  it  can  be  seen  that,  for  com¬ 
patibility  with  equation  (4),  the  quantity  oop  must  take  the  form 
o>p  =  as  +  b.  When  equations  (5)  are  substituted  into  equation  (4)  the 
results  are 


* 


9's  =  ~s^'e 


(6) 
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Elimination  of  cp*  from  the  above  equations  leads  to  the  Tricomi 
equation  • 

t'ss  +  =  0  (.7) 

Nov  equation  (4)  may  he  written  in  the  form  of  equation  (3)  by  the 
introduction  of  the  independent  variable  s*  defined  by 


are  obtained.  The  relation  between  s  and  s*  is  shown  in  figure  k 
and  tabulations  of  s,  q,  l/q,  l/p,  s*v  and  are  given  in. 

tables  I  and  II.  ■  * 
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EXACT  EQUATION  OF  STATE  LEADING  TO  THE  TRICOMI  EQUATION 


This  reduction,  which  consisted  essentially  of  changes  in  the 
dependent  and  independent  variables  and  an  approximation  of  —  ~  ^2 

p2 

to  the  second-order  term  in  its  power-series  expansion  in  terms  of  s*, 
may  he  regarded  as  producing  the  exact  equation  of  flow  for  a  gas  having 
a  somewhat  altered  pressure-density  relation.  (This  would  correspond 
to  viewing  the  Karman-Tsien  or  the  Chaplygin  approximation  for  the  holo¬ 
graph  relations  as  the  exact  equations  for  a  gas  having  a  pressure- 

density  relation  of  the  form  p  =  A  +  Using  Bernoulli’s  equation  in 

differential  form,  the  equation 


(10) 


is  obtained,  where  v  =  l/p.  If  the  approximation  m*(s*)  for  — -J— 

P2 

is  now  used,  a  differential  equation  for  v  is  obtained  whose  solution 
yields  the  relation  for  p  as  a  function  of  s*  for  the  new  gas. 
Although  the  equation  is  nonlinear,  it  is  of  the  Riccati  type,  and  can 
be  solved  either  by  expressing  v  as  a  power  series  in  s*  or  in  closed 
form  by  making  a  suitable  substitution  for  the  dependent  variable.  Since' 
the  behavior  of  the  flow  near  s*  =  0  is  desired,  the  power-series 
expansion  was  developed  and  it  was  found  that  the  specific  volume  could 
be  expressed  as 


I  =  1.5774  -  2.4883s*  +  8.6354(s*)2  -  30.175(s*)3  +  88.024(s*)4  - 

223.32(s*)5  +  .  .  .  (H) 

The  first  four  terms  agree  with  a  similar  expansion  of  l/p  for  the 
isentropic  pressure-density  relation  with  y  =  1.4. 

Using  the  relation  between  s*  and  q  (equation  (ll))  and  a 
Maclaurin  expansion  for  l/q,  the  series  expression  for  the  velocity  is 

i  =  1.0954  +  1.7280s*  +  i.7199(s*)2  -  3.934o(s*)4  +  6.7097(s*)5  - 


9.5639(s*)6  +  .  .  . 


(12) 
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From  this  and  Bernoulli's  equation  an  expression  may  he  obtained  for  the 
pressure  as  follows: 

P  -  Ps  =  0.8333s*  1.31^5(s*)2  +  0.6912(s*)3  -  0.38l6(s*)4  + 

1.0939(s*)5  -  1.39W(s*)6  +  .  .  .  (13) 

Equations  (ll)  and  (13)  then  constitute  the  pressure-density  rela¬ 
tion  leading  to  the  Tricomi  equation.  These  expressions  give  agreement 
with  the  usual  isentropic  relations  up  to  the  third  derivative  of  p 
with  respect  to  p  at  the  sonic  velocity. 

However,  the  use  of  a  power  series  to  obtain  l/p  is  hampered  by 
the  fact  that  the  series  converges  very  slowly.  As  an  alternative 
method,  the  Riccati  equation  may  be  solved  in  closed  form.  As  usual 

this  involves  a  substitution  of  the  form  v,=  i  -  —  in  equation  (10), 

P  u 

which  yields  '  : 


\ 

(l  +  cs*)^u”  =  d  X  s*u 


where  c  =  1.2304  and  d  =  9-4204  are  substituted  for  convenience. 
Letting  - 


u  =  (l  +  cs*)w(r]) 


where 

*  .  \ 

d1/^*  *;  • 

T)  =  - 

1  +  CS* 

I 

leads  to  the  differential  equation  w^  -  tjw  =  0.  This  is  the  well- 
known  Airy  equation  and  a  linear  combination  of  its  solutions  together 


* 
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with  boundary  conditions  on  l/p,  that  is,  at  s*  =  0  and  i  =  1.57744, 
provides  the  result 

1  =  1.2304 

P  1  +  1.2304s* 

2.4176  [jg/3l-2/3(l>  - 

(1  +  1.2304s*)2  _  o.1963I1/3i_1/3U)] 

where 

i  =  I.  i«6u8s»  W2 

\1  +  1. 2304s*/ 


and  Iv(l)  is  a  Bessel  function  with  imaginary  argument  in  the  notation 

of  Watson  (reference  15).  The  quantity  ^5  becomes  imaginary  for  super¬ 
sonic  values  of  speed  (s*  <  0),  and  if  |  is  defined  by 

£  =  (  -1.6118s* — \3/2  re]_a^^on  for  i/ p  in  the  supersonic  range  is 

Vl  +  1.2304s *) 


n 


1  =  1.2304  + 

P  1  +  1.2304s* 


2.4176 

(1  +  1.2304s*)2 


[l2/3j_2/3(T)  -  0.196312/3J2/3(t2) 

j^|1/3ji/3(T)  -  0.1963|1/^3J_1/3(i)_ 


where  JV(S)  is  the  Bessel  function  of  the  first  kind.  The  func¬ 
tions  q  and  M  may  also  be  computed  in  the  same  manner  and  they  are 
found  to  have  the  form 
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_ -0.1667 _ _ . 

(1  +  1.2304s*)[|1^3i1/3(|)  -  o.l96361^3l_1/3(ft)]  ' 

o  .  9.4204s*p2 

M  =  1  -  - - - 

(1  +  1.2304s*)5 

Thus  again  the  state  quantities  have  been  found  for  the  present 
hypothetical  gas  in  terms  of  s*.  This  pressure-density  relation  has 
been  compared  graphically  with  the  usual  isentropic  relation  for 
7  =  1.4  (see  fig.  5).  The  curves  are  almost  coincident  in  the  range 
p  =0.42  to  p  =  0.72.  This  corresponds  to  a  Mach  number  range  of  1.35 
to  O.85.  For  values  of  s*  near  zero,  the  curves  plotted  for  Mach  num¬ 
ber  as  a  function  of  s*  for  the  two  pressure-density  relations  agree 
quite  closely  (see  fig.  6).  For  large  values  of  s*,  however,  the 

approximation  curve  departs  from  the  usual  isentropic  curve  p  =  P^/7, 
since  the  slope  of  the  hypothetical  pressure-velocity  relation  decreases 
more  rapidly  than  that  for  the  isentropic  curve  (see  fig.  5). 


H0D0GRAFH  SOLUTIONS  YIELDING  ACCELERATED- 
DECELERATED  CHANNEL  FLOW 


Solutions  for  the  Tricomi  equation  (7)  which  produce  an  unsymmet- 
rical  accelerated-decelerated  channel  have  been  given  by  Tomotika  and 
Tamada  (reference  3)>  Falkovich  (reference  1 6),  and  Ehlers  (refer¬ 
ence  13).  Consider  a  solution  which  has  the  form 


to(B,0)  =  ijj-  (^e2  +  |  s3  _  ej1^3  -  (^e2  +  ^  s3  +  e)1^3  (i4) 


Tomotika  and  Tamada  noticed  that  this  function  remained  a  solution 
when  (0  -  ik)  is  substituted  for  0.  In  fact,  if  the  function 

-  >/3 

bj  _  1 fl  _  in  1  r 

9 


♦i(s ,6)  =  Re  -i- 
"L 


fie  -  i\)2  +  |  e3  .  (e  -  ix)J 


1 1  ^ 

-  iA,)2  +  i  s3  +  (e  -  ix ) 


(15) 


< 
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is  considered,  a  stream  function  results  which  has  a  branch-point  sin¬ 
gularity  of  the  one-half  power  at  s  =  ^  =  sQ  and  0  =  0.  This  « 

branch  point  in  the  subsonic  stream  may  be  said  to  " generate"  a  channel 
flow  with  locally  supersonic  regions  having  a  Mach  number  at  the  nar¬ 
rowest  cross  section  determined  by  A,.  Of  the  other  properties  of  the 
function  it  may  be  observed  that,  when  s  >  sQ  and  0=0,  then 
^  =  0.  It  would  appear  that  \|f]_(s,@)  has  a  line  of  one- third- order 
branch  points  along  s  =  0.  However,  this  branch  point  may  be  elimi¬ 
nated  by  writing  the  function  in  the  form 


3 
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where 


Pi  =  [|a  cos  |  -  «)2  +  (a  sin  |  +  i 

a.^^s3-  x2  *  e2)2  +  ia2e2 


ft  =  tan-1 


a  sin  ®.+  X 


a  cos  |  -  0 


it  s3  -  x2  +  e: 
9 


Since  does  n0,t  vanish  (the  choice  of  sign  for  the  one-half -order 

"branch  point  is  such  that  this  is  insured)  it  can  he  seen  that  the  sin¬ 
gularities  along  s  =  0  are  removable.  The  observation  may  be  made 
that 


ix(s,e)  +  t1(s,-e)  =  f(s,e 


is  a  symmetrical  function  of  0,  and  it  is  this  function  which  shall  be 
used  as  the  stream  function.  The  function  $  has  "been  plotted  in  fig¬ 
ure  7  for  X  =  0.050912.  This  corresponds  to  s  =  0.l8  at  the  center 
of  the  channel  or  a  maximum  Mach  number  on  the  center  streamline  of  0.86, 

For  sufficiently  large  negative  values  of  s  ( supersonic  flow) ,  ||- 

becomes  zero  at  a  point  above  the  0  =0  axis  as  well  as  at  0  =  0. 

This  indicates  the  flow  in  the  supersonic  region  will  have  streamlines  / 
along  which  the  velocity  gradient  changes  sign  more  than  once.  This 
can  be  shown  from  the  following  analysis.  A  series  expansion  of 

for  constant  s  about  0=0  has  the  form  ! 


*e  -  (*ea) 


0  +  0 
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From  the  holograph,  it  is  apparent  that  for  s  >  0.l8  the  derivative 
$0  is  positive  for  0  >  0.  Along  a  streamline  for  which  does 

not  vanish  for  a  value  of  9  >  0,  the  value  of  then  will  be 

positive  in  the  neighborhood  of  the  line  0=0  since  0  for 

®=0.  The  limiting  streamline  of  those  along  which  the  velocity  gra¬ 
dient  does  not  change  sign  more  than  once  is  the  streamline  which  passes 
through  the  point  on  the  s-axis  for  which  fee  vanishes.  This  was 

found  to  be  f  =  0.762  for  s  =  -0.132.  For  s  <  -0.132  and  hence 
$  >  0.762  the  value  of  for  9  >  0  in  the  neighborhood  of  the 

line  0=0  is  negative;  for  s  >  -0.132  and  $  <  0.762,  is 

positive.  (See  fig.  8(a).)  The  resulting  effects  on  the  constant- 
velocity  curves  in  the  physical  plane  are  shown  in  figure  8(b).  It  is 

apparent  from  the  figure  that  ||  =  0  does  not  occur  along  the  charac¬ 
teristics  passing  through  s  =  0  and  0=0.  This  remark  is  made  at 
this  point  since  in  the  limiting  case  as  \  — »  0  this  characteristic 
does  become  the  locus  of  the  changes  in  sign  of  velocity  gradient.  This 
is  discussed  in  more  detail  in  the  section  Limiting  Case  When  Sonic 
Velocity  Is  Just  Reached  on  Center  of  Streamline. 

The  breakdown  of  isentropic,  potential  flows  is  usually  associated 
with  the  appearance  of  limit  lines.  Some  remarks  have  already  been  made 
about  criterions  for  breakdown,  but  here  it  will  be  sufficient  to  remark 
that  Tollmein  (reference  17)  and  Tsien  (reference  18)  have  shown  that 
inviscid  potential  flows  cannot  be  continued  across  limit  lines.  Limit 
lines  occur  when  the  Jacobian  of  the  transformation  from  the  holograph 
to  the  physical  plane 


(20) 


becomes  zero.  This  may  only  occur  for  s  <  0  since  s,  and  $s 

are  real  and  the  Jacobian  is  always  positive  for  s  >  0.  Geometrically, 
a  necessary  and  sufficient  condition  that  limit  lines  occur  is  that  the 
streamlines  and  the  characteristics  be  tangent  in  the  hodograph  plane, 
provided  the  derivatives  and  $s  remain  finite.  As  a  matter  of 

fact  in  this  case  the  limit  lines  are  the  locus  of  the  points  of  tangency. 

The  hodograph  of  the  function  f  was  tested  graphically  for  limit 
lines  by  seeking  points  of  tangency  of  the  streamlines  using  a  "portable" 
characteristic  curve,  drawn  on  tracing  paper  and  moved  across  the  super¬ 
sonic  field.  From  figure  7  it  is  apparent  that  no  limit  lines  occur  for 
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that  portion  of  the  s,0  field  which  has  "been  computed.  Using  this 
test,  it  is  not  certain  that  limit  lines  do  not  appear  somewhere  else 
in  the  supersonic  portion  of  the  flow,  but  since  the  choice  of  bounding 
streamline  has  so  far  been  left  open  this  is  not  an  important  problem. 


CALCULATION  OF  PHYSICAL  COORDINATES 


To  obtain  the  coordinates  of  the  flow  in  the  physical  plane  the 
condition  of  irrotationality  and  continuity  may  be  put  into  the 
following  form: 


dcp  =  u  dx  +  v  dy 

-  dt  =  (-v  dx  +  u  dy)i 
P 

These  equations  may  then  be  written  as 


Expressing  cp  and  \|r  as  functions  of  s  and  0,  and  then  replacing  ' 
the  partial  derivatives  of  cp  with  those  of  \| r  by  equation  (4),  the 
following  differential  expression  for  the  complex  position  coordinate 
is  obtained: 

+  i  ts)  ds  (22) 

With  this  formula  it  is  possible  to  show  that  if  t  is  an  even 
function  of  0  the  flow  pattern  in  the  physical  plane  is  symmetric. 

For  this  purpose,  it  is  first  of  all  necessary  to  prove  that  the  line 
0=0  in  the  hodograph  plane  maps  into  a  straight  line  in  the  physical 
plane.  Setting  0=0  and  d0  =  0  in  equation  (22)  leads  to 


dz  =  -  eie 

q 


(-"A  +  jj  *e) 


d  e  + 


* 


dz  =  i(sa)n2^ 

q\ 


1  Ye 


\=o 


ds 


(23) 
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The  stream  function  $  in  equation  ( 19)  is  an  even  function  of  0}  hence 
$0=0  for  0  =  0.  Then 

d2 '  S  **L ds  -S  d*L  <2« 

This  is  purely  imaginary;  and  0=0  then  is  seen  to  map  into  the 
y-axis.  If  the  integration  for  the  streamlines  is  carried  out  along 
each  streamline  from  the  point  0=0,  the  obvious  differential  relations 


dty  =  0  =  \|fg  ds  +  d0 


lead  to  the  simple  formula,  where  \|r  =  Constant, 


(25) 


Since  (see  preceding  section)  the  Jacobian  of  the  transformation  from 
the  s,0-  to  the  x,y-plane  is  never  zero,  and  since  it  may  be  represented 
as 


* 


Equation  (25)  becomes 


\ 


v 


^Constant 
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x 


pq  cos  0  J 


d  0 


t=Constant 


del 

/  ty=Constant 

Nov  for  f  In  equation  (19),  the  relation  p(e)q(e)  =  p(-0)q(-0)  along 
any  single  streamline,  since  $s  and  f  are  even  in  0,  implies  that  the 
same  value  of  s  occurs  for  0  and  -0.  Thus  for  x  the  integrand  is 
a  product  of  four  even  functions  of  0:  J,  pq,  $s,  and  ^  cos  0.  There¬ 
fore,  x  itself  is  odd.  In  a  similar  manner,  y  is  seen  to  he  an  even 
function  of  0.  Since  x  is  an  odd  function  of  0  for  a  given  value 
of  s  and  y  is  an  even  function  of  0,  then  the  channel  is  symmetric 
about  the  line  0=0. 

(  \ ' 

The  mapping  from  the  hodograph  plane  back  to  the  physical  plane 
may  be  accomplished  by  integrating  equations  (2l)  or  (22)  along  either 
constant  '  \|r  or  constant  s  lines.  Since  either  one  of  these  inte- . 
grations  must  be  carried  out  by  numerical  or  graphical  methods,  the 
choice  depends  on  the  method  involving  the  least  amount  of  labor. 


y  =  - 


P0 


pq  sin  0  J 


If  the  integration  is  performed  along  a  streamline,  then  from  equa¬ 
tion  (21)  the  physical  coordinates  are  found  by  integrating 


dz  =  i  ei0  dcp 

q 

If  cp  were  known  as  a  function  of  s  and  0  along  each  streamline 
this  would  yield  a  very  simple  method  for  computing  the  flow.  However, 
the  differential  d9  may  be  replaced  by  the.  differential  of  <p'  with 
an  additional  term  involving  i r  and  o)p.  Multiplying  the  second  of 
equations  (5)  by  ds  and  the  third  by  d0  and  adding  the  two  relations 
lead  to 

1  to-, 

dcp*  = - dcp  +  -  t  d0 

ds 
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Thus  the  physical  coordinates  are  found  by  integrating 


dz  = 


<1 


+  —  eie\|r  d0 

q 


(26) 


along  \(r  =  Constant. 

Nov  the  function  qp*  can  be  computed  from  equations  (6)  when  \|r' 
is  known.  Differentiating  ti'  with  respect  to  s  gives 


n - - — ; — -  i  -2/3 

_  O 

Re  m 

i|j(0  -  iX)2  +  i  S3  -  (0  -  iX)j 

2s^ 

3 J(0  -  iX)2  +  ^  s3 

L 

L 1  9 

r*  .  _ _  -i 

-2/3 

—  “1 

|J(0  -  iX)2  +  £  s3  +  (0  -  iX) 

2s2 

> 

12  y  J 

3J0  -  iX)2  +  J  s3 

L<  9  J 

This  expression  may  be  simplified  by  multiplying  the  first  term  inside 
the  braces  by  unity  in  the  form 


“I  2/3 
+  (0  -  iX) 


12/3 

+  (0  -  iX) 


and  by  multiplying  the  second  term  by  a  similar  factor  with  the  posi¬ 
tive  sign  before  ( 6  -  iA,)  changed  to  minus.  This  leads  to 
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The  arbitrary  function  of  s  results  from  the  integration  with  respect 
to  0.  That  f(s)  =0  will  be  readily  apparent  for  the  functions  i|r ' 
and  -<p'  must  obey  the  relation 


«PS'  =  -^0 


2k 
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Thus,  differentiating  equation  (28)  with  respect  to  s  yields 


Employing  the  same  procedure  as  on  >|rsf  leads  to 


The  expression  on  the  right  side  can  be  readily  recognized  as  -s^1  as 
required.  The  potential  <p*  corresponding  to  given  by  equation  (17) 

is  then 
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If  for  constant  s  the  function ..  is  an  even  function  of  0 , 
then  qf  is  an  odd  function.  This  is  apparent  from  equation  (6)  above, 
since  differentiation  of-  an  odd  function  of  G  with  respect  to  9 
results  in  an  even  function.  Thus,  the  conjugate  potential  corre¬ 

sponding  to  $  is 

«’■  =  q>’(s,e)  -  cp' (s,-e) 


Here  use  has  been  made  of  the  fact  that  0^  an  even  function  of  G . 
In  order  to  perform  the  actual  integration,  the  streamline  and  con¬ 
stant  cp1  lines  must  be  plotted  in  the  hodograph  plane.  Then  for  each 
streamline  the  quantities  (o^i/q)  cos  G  and  (o)q/q)  sin  6  are  plotted 
against  and  also  against  0.  The  integrations  then  are  performed 

with  the  aid  /of  a  planimeter. 

A  symmetric  channel  flow  has  been  computed  for  sQ  =  0.l8 
(\  =  0.0509117);  which  corresponds  to  a  Mach  number  of  0.86  at  the 
center  of  the  throat  section  (see  figs.  8(a)  and  8(b)). 


VELOCITY  DISTRIBUTION  ON  CENTER  STREAMLINE  f  =  0 


The  streamline  $  =  0  must  be  treated  separately.  The  position 
coordinates  are  computed  by  integrating  equation  (22),  which  for  this 
case  reduces  simply  to 


dx 

ds 


N 


2  6 
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since  $s'  =  0  and  d0  5  0.  The  derivative  with  respect  to  0  is  found 
by  differentiating  the  stream  function  in  the  form  of  equation  (15). 

This  leads  to 


r10 


Re 

3^ 


[j(0  -  iX)2  +  i  S3  -  (0  -  iX)] 


I1/ 3 


^(0  -  ix)2  +  i  s3 


or 


jj(0  -  iX)2  +  ^  s3  +  (0  -  iX)j 
^(0  -  iX)2  +  ^  s3 


(31) 


"  -  3^ 


Q(0  -  iX)2  +  |  s3  -  (0  -  iX)]l/3 


(0  -  iX)2  +  |  s3 


$1/3' 


^(0  -  iX)2  +  i  s3  j( 0  -  iX)2  +  ^  s3  -  (0  -  ix) 


1/3 


(32) 


Finally, 


+10  =  " 


1 

30)30- 


ny/3s 


(33) 
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With  this  expression  substituted  for  the  physical  coordinates 

are  found  by  integrating 


dx 


-  r3(i) 


4/3  s3/2a>,_ 


cos 


qao 


(S-f) 


ds 


(35) 


h  \2/3 

Unfortunately,  a0  =  0  at  s  =  If  M  and  the  quantity  on  the  right- 
hand  side  becomes  infinite  so  that  numerical  integration  in  this  form  is 


impossible.  This  can  be  remedied  by  noting  that 

ds  ~~ 

this  substitution  for  Oq,  equation  (35)  yields 


With 


2  sa. 


dx 

d& 


=  -2 


Vb  s5A 


\q 


cos 


(i-t) 


(36) 


*  \ 


The  integration  can  be  readily  performed  by  plotting  the  right-hand 
side  against  6  and  using  a  planimeter.  The  velocity  distribution 
along  $  =  0  as  well  as  along  the  other  streamlines  is  plotted  in  fig¬ 
ure  9  for  X  =  0.050912. 


CURVES  OF  CONSTANT  VELOCITY 


The  curves  q  =  Constant  in  the  physical  plane  are  found  by  inte¬ 
grating  the  relation 


dz 

d6 


^2fs 


+  —  % 


pq 


q=Constant 


(37) 


To  study  the  nature  of  the  curve  in  the  neighborhood  of  the  vertical 
axis  of  symmetry  x  =  0  (0=0  and  s  <  0.l8)  it  is  convenient  to  use 

a  power-series  expansion  about  0=0.  The  coefficients  of  such  a  series 
are  readily  obtained  by  differentiating  successively  the  expression  on 
the  right  with  respect  to  0  and  then  setting  0=0.  The  expansion 
for  z  becomes  . 
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which  for  0  =  0  reduces  to 

K  =  2c  =  (  q$ee  \ 

A2  \p%^s2/e=0 


The  sign  of  the  curvature  at  the  point  0=0  depends  upon  the  sign 
of  C  which  in  turn  depends  on  the  sign  of  $ 00*  Now,  the  second 

derivative  becomes  zero  for  s  =  -0.132  and  0=0.  For 

s  >  -0.132,  the  quantity  t00  is  greater  than  zero  and  the  curvature 
of  the  s  =  Constant  lines  is  convex  to  the  center  streamline  of 
the  channel  (see  fig.  8(b)).  For  s<  -0.132,  the  quantity  is 

less  than  zero  for  0=0  and  the  curvature  of  the  constant  s  line 
is  concave  to  the  center  streamline  of  the  channel.  .  The  curves 
s  =  Constant  for  0.l8  >  s  >  -0.132  show  a  minimum  along  x  =  0, 
while  the  curves  s  =  Constant  for  s  <  -0.132  exhibit  a  maximum  at 
x  =  0,  and  two  minimums  symmetrically  placed  about  the  y-axis.  Along 
the  streamlines  which  cut  these  latter  curves,  the  fluid  is  accelerated 
to  a  maximum  velocity  and  then  is  decelerated  slightly  to  a  minimum 
velocity  at  the  throat.  The  fluid  is  accelerated  to  the  maximum  velocity 
before  being  decelerated  finally  to  subsonic  velocities  again.  This 
behavior  is  apparent  from  the  velocity  distribution  along  the  streamlines 
$  =  0.832  and  0.920  in  figure  9- 

The  curves  s  =  Constant  for  s  >  0. 18  which  intersect  the  x-axis 
at  90°  are  concave  to  the  incoming  flow  in  the  neighborhood  of  the 
straight  streamline  at  the  entrance  side  of  the  vertical  axis  of  sym¬ 
metry  and  convex  to  the  incoming  flow  on  the  exit  side.  Further  from 
the  center  streamline  and  toward  the  wall,  however,  the  curvature  changes 
and  the  constant-velocity  lines  become  convex  to  the  incoming  flow  on 
the  entrance  side  and  concave  to  the  incoming  flow  on  the  downstream 
side  of  the  vertical  axis  of  symmetry  (see  fig.  8(b)). 

Since  the  stream  function  has  a  branch  point  at  s  =  0.l8  and 
0=0,  the  above  analysis  cannot  be  applied  to  the  curve  of  constant 
velocity  s  =  0.l8.  To  investigate  the  nature  of  this  curve  in  the 
neighborhood  of  the  center  streamline,  it  is  necessary  to  obtain  the 
expansion  of  equation  (37)  for  small  values  of  0.  The  derivative  tie 
is  found  from  equation  (33).  Using  the  fact  that  ft(-0)  =  it  -  fi(0) 
and  cd(-0)  =  2it  -  co(0)  leads  to 
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J L 

3%a 


-  Pt  cos  /----•£■ 
a  1  \3  2  6, 


\3  2  6y 


where  Si,  a,  and  cd  are  defined  in  equations  (l8).  In  a  similar  way 
by  using  equation  (27)  for  'l'ls'(0)  the  following  equations  are 
obtained: 


*S  =  StEls‘(0)  +  +ls'(-0O  -  k 


Letting  s  =  0.l8  in  equations  (39)  and  (4o)  leads  to  some  simplifica¬ 
tion  since  for  this  case  ^  s^  =  X 2.  Accordingly,  for  values  of  9 
positive  and  near  zero, 

a  -  e1^  ♦  kx^k  *  (2 iff)  ^  ♦  o(e5/2)  1 


=  tan* 


(fH 


1  ^ 

SI  =  tan-1  *- ...  ■ 

.  a,2  +  ©2 


f*(*J 


|1/2  +  oW‘) 


.2  ...  o2 


6i  ■  e)  +  x 
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The  approximate  values  of  the  above  quantities  lead  finally  to 


■Jr =  A 

6  n|3 

=  0. 729020"  +  ; 

5.7280  +  o(e3/2)  +  .  .  !] 

(42) 

u.-^i 

'4\1/3  0-1/2 

2a  ri/6  +  _l] 

0l/2 

*s  2  ' 

<9/ 

36  \9J  a>j_X 

=  -0.30930' 

■1/2  _  5.87701/2  +  o(03/2) 

(43) 

With  these  expansions  and  equation  (37),  the  slope  of  the  curve  s  =  0.l8  m 

in  the  neighborhood  of  0=0  is  given  by 

&  =  1.948  -  39.530  +  o(e2) 

dx 

The  numerical  value  of  the  slope  at  0=0  is  finally 

^  =  1.948 


Now  the  curvature  is  proportional  to  the  second  derivative 
for  the  present  case  is  given  by 


d2y 

— which 
dx2 


2£  «  e1/2  ♦  o(e) 

d 0\dx/  dx 


The  curvature  of  the  s  =  0.18  line  is  seen  to  be  zero  where  it  cuts 
the  center  streamline. 
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LIMITING  CASE  WHEN  SONIC  VELOCITY  IS  JUST  REACHED 
ON  CENTER  OF  STREAMLINE 


Setting  X  =  0  in  the  stream  function  in  equation  (15)  leads 
directly  to  the  stream  function  considered  for  the  flow  through  a  Laval 
nozzle  in  reference  13.  From  this  solution  it  is  possible,  however,  to 
construct  a  channel  having  locally  supersonic  regions  hy  the  appropriate 
choice  of  tranches.  Tomotika  and  Tamada  (reference  3)  constructed  a 
channel  having  M  =  1  at  the  center  streamline  using  a  solution  in 
the  \|r,cp-plane.  Their  solution  for  this  limiting  case  can  be  shown  to 
be  the  same  as  equation  (l4)  with  X  =  0  when  expressed  in  the  hodo- 
graph  variables.  Now  equation  (6)  when  expressed  in  cp*  and  V  as 
the  independent  variables  becomes 
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of  s  and  0.  Setting  the  factor  z  -  2\x  equal  to  zero  and  elimi¬ 
nating  cp*  and  z  with  the  aid  of  equation  (45)  leads  to 


(■ r )3  + 


0 


The  resulting  function  t  is 


Similarly,  the  relation  z  +  p  =  0  yields 


(46) 


(Mr1)3 


3@ 

8 


0 


whose  solution  leads  to  the  stream  function 


Equations  (46)  and  (47)  express  the  two  tranches  of  the  stream  function 
which  go  to  make  up  the  accelerated-decelerated  flow  through  a  nozzle. 
Equation  (46)  may  he  obtained  from  equation  (47)  by  inserting  unity  in 
the  form  of  e2jti  inside  of  the  parentheses  on  the  right-hand  side  of 
equation  (47)  and  taking  the  real  part.  For  the  flow  computed  by 
Tomotika  and  Tamada,  figure  10  shows  the  regions  of  the  hodograph  plane 
in  which  the  two  branches  represented  by  equations  (46)  and  (47)  were 
used. 


To  obtain  a  flow  which  is  symmetric  about  the  line  0=0  as  well 
as  about  the  straight  streamline  the  stream  function  given  by 
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f '(e)  =  ffHe)  +  t(-0)| 


m 


shall  be  considered.  In  this  formula  \|r(0)  is  given  by  equation  (47) 
and  i|r( -0)  "by  equation  (46)  if  0  >  0  and  02  >  ^  (-s)3.  in  the  region 

2  3/2 

between  the  characteristic  0  =  (-s)  ,  \|/(0)  and  \|r(-0V  are  both 

given  by  equation  (47).  The  choice  of  the  factor  2/3  was  made  so  , 
that,  in  the  region  of  the  upper  quarter  plane  0  >  0  and  s  >  0  and 

?  3/2 

the  section  in  the  left  quarter  plane  where  0  >  £  (- s)  ,  the  stream 

function  is  simply  given  by 


?  =  — 
o>i 


(«  -  >F^P)1/3  *  («  +  ^  7  ^  73) 


1/3 


(49) 


The  physical  plane  for  this  portion  o^  the  hodograph  has  already  "been 
computed  by  Ehlers  (reference  13) • 

The  relation  for  the  stream  function  in  the  region  between  the 
2  3/2 

characteristics  6  =  ±j(-s)J/  is 


,  >l/2|  5(0)  8(-0) 

(-s)  COS  -1 — -  +  cos  — —  ' 

1  3  3 


3 


(50) 


where 


5  =  cos“^ 


(-B) 


3/2 


5(0)  =  n  -  8(-0) 


0  >0 
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Hence,  with  the  aid  of  a  simple  trigonometric  identity  the  stream 
function 


27/33l/6 

3u>i 


COS 


(51) 


is  obtained.  The  streamlines  in  the  holograph  plane  for  the  entire 
flow  given  hy  equation  (k8)  are  shown  in  figure  11(a). 

To  compute  the  coordinates  for  the  curves  of  constant  speed  near 
the  y-axis  the  relation 


was  employed.  By  a  MacLaurin  expansion,  for  small  values  of  0,  the 
x  coordinate  can  he  computed  from 


x 


+  0)]_ 


fQe  +  o{e3) 


where  is  the  value  of  the  stream  function  when  0=0.  The  y 

coordinate  is  computed  from 


For  small  values  of  0,  this  is  approximated  hy 


y  -  =  —(f  -  fo)  +  °(0cp) 


where  cp  =  0  at  0=0.  The  flow  in  the  physical  plane  is  shown  in 
figure  11(h). 
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The  following  paragraphs  will  give  a  discussion  of  the  flow  pattern. 
The  derivatives  $0  and  $s  of  the  stream  function  in  equation  (51) 

o  3/2 

become  infinite  at  the  line  of  branch  points  6  ±~  (-s)  =  0.  These 

curves  are  the  characteristics  in  the  hodograph  plane  which  are  cusped 
at  the  origin  of  coordinates.  As  these  characteristics  are  approached 
along  the  branch  of  the  streaip  function  given  by  equation  (49)  tyg 
and  remain  finite;  thus  there  is  a  jump  in  the  first  derivatives 

at  the  characteristics  6  ±|  (-s)^/2  =  0.  Furthermore  along  the 

f  =  0  streamline,  $s  and  have  the  following  values: 


= 


(2/3) 


1/3 


CD]_S 


The  quantity  is  seen  to  be  singular  at  the  sonic  point. 

The  effect  of  the  above  singularities  on  the  character  of  the  flow 
will  now  be  investigated.  From  the  equation  transforming  the  solution 
from  the  hodograph  to  the  physical  plane,  namely, ■ 


=Constant 


ds 


and  the  relation 

p  oo;i2P 

ds  =  cdt  ds*  =  -  -  dq 

-1-  q 

the  velocity  gradient  along  a  streamline  can  be  expressed  in  the 
following  form: 


/d q \  =  _  q2  /  f£  \ 

W_4-k  a^P\$s2  +  sf02/ 


m12p2J 


% 
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where  J  is  the  Jacobian  of  the  transformation  from  the  hodograph 
plane  to  the  physical  plane.  As  the  Mach  line  corresponding  to  the 

characteristic  0  =  |  (-s)3/2  is  approached  from  the  sonic  point  of 

the  streamline,  the  velocity  gradient  approaches  a  finite  value  since 
and  J  are  both  finite.  As  the  same  Mach  line  is  approached  from 
the  direction  of  the  0  =  0  line,  the  quantities  J  — ►  00  and 

i| iQ  — ■>  00  in  such  a  way  that  ^  is  finite  but  differs  in  value 

from  when  approached  from  the  other  direction.  The  line 

\dZ/$= k 

J  =  00  then  corresponds  to  a  discontinuity  in  the  velocity  gradient 
along  the  streamline.  For  the  center  streamline  $  =  0,  the  velocity 
gradient  becomes 


.  n  2 

_ _ 5 -  (52) 

<ix  CD]_2ps$g 


Now  for  f  =  0,  °c  i  as  s  — *  0.  Since  s$0  has  opposite  signs 

on  the  two  sides  of  the  slit  corresponding  to  $  =  0  in  the  hodograph 
plane,  and  approaches  a  finite  nonzero  value  as  s  — >  0,  there  is  a 
jump  in  the  velocity  gradient  given  by 


/dq\  =  (2q*)2 

W"°  '  <1.4532)V(f)l/3 

where  q*  and  p*  are  the  critical  velocity  and  density,  respectively. 
To  study  the  curvature  of  the  streamlines,  the  relation 

(a)  . .  ai  Mt.  * 

\d®/^=Constant  ^  V  ^s  / 

is  considered.  In  a  manner  similar  to  the  computations  of  the  velocity 
gradient  along  the  streamline,  the  following  result  is  obtained  for  the 
curvature : 


(53) 
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This  expression  is  discontinuous  at  the  line  J  =  «>,  since 


is 

j 


approaches 


different  values  as  the  characteristic  (j  =  »)  is  approached  from  the 
two  directions  along  the  streamline.  Thus,  the  Mach  line  (j  =  oa)  is  a 
curve  at  which  the  curvature  of  the  streamlines  as  well  as  the  velocity 
gradient  is  discontinuous. 


SOME  OTHER  SOLUTIONS  OF  THE  TRICOMI  EQUATION 


Tricomi  (reference  11)  has  shown  that  a  single-parameter  family  of 
solutions  of  the  differential  equation 


+  sf*00  =  0 


(54) 


is  given  "by 


r 


c^F 


(■“' 


n  + 


1  2 
3j  3j 


— 1 
2a  J 


(55) 


where  a 


Now  it  is  possible  to  derive  from  this  formula 


some  solutions  in  closed  form.  In  place  of  the  hypergeometric  function 
above,  consider  the  following  related  functions  (see  p.  286,  Whitaker 
and  Watson  (reference  19)): 


t’i 


1 


n 


=  a11 


(-  ^r/3 


(n 


+ 


— ,2n  + 

3 


4.  _2a_\ 
3*  a  +  ej 


and 


(56) 


In  \|r' ,  let  n  +  — 

3 

solutions  "become 


m;  and  in 


let  -n  +  — 
.3 


m.  1  Then  the 
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♦  •  =  cP~2M - 1/,3  F(m  -  m;  2m;  — 

l,m  \  a  +  6/  \  3  a  +  6/ 


*'2,m  =  a ' 


■m+ 


l/3(-  TTsf  1/3  F("  '  5;  “J  2n;  77%) 


(58) 


Thus  it  is  apparent  that 


V .  =  tr2111"1  if'  , 

v  l,m  y  2,m* 


(59) 


for  a  given  value  of  m. 

In  order  to  see  how  these  solutions  may  be  summed,  consider  the 
hypergeometric  function  written  in  series  form;  namely. 


r(m)r| 

(*-l) 

00 

_  I*  (2m) 

c  - 

r(m) 

z_ 

n=0 

T(2m) 

d 

r(m) 

dzm_1 

_  r(2m) 

d 

r  +  m  -  ijr(m  +  n)zn 
r(2m  +  n)n! 


n=0 


r  +  m  -  ^zn(m  +  n  -  l)(m  +  n  -  2)  .  .  •  (n  +  l) 


("  •  t)( 


2m  +  n  -  l)  1 


(n  +  m  -  i) 


,n+m- 1 


r(m)  dz111-1 


Pz 


rff-  a  -  d^-i 

L  t-i) 

v - v1 -  n-0 


2m- 1 
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4-1 


The  integrations  involved  in  the  above  formula  can  he  very  easily  per¬ 
formed.  Putting  m  =  1  in  equation  (60)  yields 


t'n  =  |[(e  +  a)1/3  -  (e  -  a)1^ 
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Putting  m  =  1  into  equation  (6l)  yields 

*  “>1/3  -  -  «>1/31 
Putting  m  =  2  into  equation  (60)  yields 

Vlf2  =  fa  +  ^)4/3--  (0  -  cr)4/^]  -  a[e  +  a)l/3  -  (6  -  0 y)l/3\ 

Now  (0  +  a)1/3  and  (0  -  cr)^3  are  separate  solutions  of  the  dif¬ 
ferential  equation.  The  functions  ty'2,1  can  constructed  from 

solutions  of  the  form  by  differentiating  with  respect  to  0. 

Such  a  process  of  successive  differentiation  yields  the  following  solu¬ 
tions  of  the  Tricomi  equation: 


'if'  1  =  ( c  +  0) 


1/3 


f  2  = 


♦  e)1^ 


(a  +  e) 


(a  +  e) 


1/3 


a2  * 


3  (a  +  e)4^3 

4 


r,-Uia  +  e) 


+  e)1^3  9  (o  +  e)k/3  (a+0)1/3 

*  7 


^  (62) 


7 


(a  +  0)1/3  -H(a+  6)V3  +  m  4  (a  +  0)T/3  - 


2  <P 


(a  +  0) 


128  ff6 


(a  +  0) 


Sl  i  <„  ♦  e)10/3 

56  0f 
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These  five  fundamental  solutions,  together  with  the  set  with  -0 
replacing  0  may  he  combined  to  form  the  solutions  corresponding  to 
m  =  1  in  equation  (60)  and  m  =  1,  2,  3;  and  4  in  equation  (6l), 

respectively.  By  equation  (59) ,  it  is  seen  that  cnjr'  cA^, 

and  are  also  solutions  of  the  differential  equation.  This  may 

he  verified  hy  substituting  them  into  the  differential  equation. 


The  functions  and  have  the  property  that  the 

point  s  =  0  and  0=0  is  mapped  into  the  point  at  »  in  the 
physical  plane.  In  fact  the  solutions  all  become  indeterminate  at  that 
point.  Thus,  these  solutions  may  be  used  to  construct  flows  having  a 
uniform  sonic  stream  at  infinity.  If  0  is  replaced  by  0  -  iX  and 


g  b3/2  = 


the  real  part  is  considered,  then  the  point  0=0  and  =  s 

into  infinity  in  the  physical  plane.  -Tomotika  and  Tama da  used  the 
stream  function 


X  maps 


=  Re  [^(0  -  iX)  -  tx(-0  +  iX)  +  i|r2(0  -  iX)  +  i|r2(-0  +  iX)] 


to  construct  the  flow  over  a  cusped  body  having  locally  supersonic 
regions. 

The  solutions  -  ty2(e>s)],  -  1^3(0, s)J  , 

-  ^(0,s)j,  and  -0,s)  -  *5(0, s)j  may  be  used  to 

construct  channel  flows  having  locally  supersonic  regions.  Let 


=  |j(0  -  ix)2  +  i  s3  .  (e  .  ix)j 


Since 


||(0  -  iX)2  +  ^  s3  +  (0  -  iX)j  n  = 


^ne3n 


the  appropriate  solutions  may  be  written 
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where  a(\)  =  ^(6  -  iX)2  +  ^  s^.  Expressing  the  real  parts  in  terms  of 
the  magnitudes  and  arguments  of  the  complex  quantities  yields 
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If  6  goes  to  zero  through  positive  or  negative  values  for 

^2/3 


S  > 


then  03  equals  2jt  or  zero,  respectively.  Hence 


{§)i,6*i/2 


and  the  stream  functions  vanish  identically. 


however , 


03 


jt  at  0  =  0  and 


For 


and  the  stream  functions  reduce  to  non-trivial  functions  of  s.  Thus, 
the  stream  functions  in  equations  (63),  (64),  and  (60)  will  yield 
channel  flows  with  locally  supersonic  regions.  It  should  be  pointed 
out  that  the  solutions  given  by  equation  (60)  for  all  integral  values 
of  m  >  1  seem  to  be  reducible  to  the  form 


V  -  02"-1  ^  *  e)1^ 

69 


The  solutions  for  higher  values  of  m  than  those  given  in  equation  (62) 
have  a  similar  form  and  may  be  combined  to  give  channels  having  locally 
supersonic  regions.  It  is  then  apparent  that  it  is  possible  to  approxi¬ 
mate  quite  closely  the  flow  through  a  given  nozzle  by  taking  a  sum  of 
these  solutions,  providing  the  Jacobian  of  the  transformation  to  the 
physical  plane  does  not  vanish. 

The  stream  functions  given  in  equations  (63);  (64),  and  (65)  with 
\  =  0  cannot  be  continued  into  the  supersonic  region,  since,  at  0=0 
and  s  =  0,  ^  and  are  both  zero  and  hence  the  Jacobian  vanishes 

at  this  point.  These  solutions,  however,  may  be  combined  with  the 
solution  in  equation  (l4). 


Brown  University 

Providence,  R.  I.,  March  13,  1950 
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TABLE  I.-  TABULATION  OF  e,  l/q,  q,  AND  l/p 


s 

i/q 

q 

l/p 

0.34 

1.45156 

1.27356 

.2  6 

1. 35002 

1.33361 

.18 

1.26118 

1. 39801 

.16  • 

1.2I065 

.80603 

1. 41545 

.14 

1.22071 

.81920 

1.43333 

.12 

1.20135 

.83240 

1.45175 

.10 

1.18252 

.84565 

1. 47078 

.08 

1.16420 

.85896 

1. 49050 

.06 

1.14636 

.87233 

1.51096 

.04 

1.12897 

.88576 

1.53223 

.02 

1.11200 

.89928 

1.55437 

.00 

1.09544 

.91288 

1.57744 

-.02 

1.07927 

•92655 

1. 60151 

-.04 

1.06345 

.94034 

1.62666 

-.06 

1.04797  • 

.95423 

1.65293 

-.08 

1.03282 

.96822 

1.68042 

-.10 

1.01796 

.98236 

1.70918 

-.12 

1.00339 

.99662 

1. 73929 

-.14 

.98908 

1.01104 

1.77084 

-.16 

.97502 

1.02562 

1.80387 

-.18 

.96120 

1.04037 

1.83848 

-.22 

.93420 

1.07044 

1.91318 

-•27 

.90152 

.  1.10924 

2.01693 

-.32 

.86986 

1. 14961 

2.13391 
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TABLE  II.-  TABULATION  OF  s, 


AND  a>L 


s 

s* 

<*1 

0.34 

0.20075 

1.16534 

.26 

.14508 

1.23307 

.18 

.  09521 

1.30080 

.16 

.08355 

1.31774 

.15 

.07782 

1. 32620 

.14 

.07217 

1.33467 

.12 

.06109 

1.35160 

.10 

.05028 

1. 36854 

.08 

.03973 

1.38547 

.  06 

.02944 

1. 40240 

.04 

.01939 

1.41933 

.02 

.00958 

1. 43627 

.00 

.0000 

1. 45320 

-.02 

7.00936 

1.47013 

-.04 

-.01851 

l.ii8707 

-.06 

-.02745 

1.50400 

-.08 

-.03619 

1.52093 

-.10 

_ . 04474 

1.53786 

-.12 

-.05311 

1. 55^*80 

-.14 

-.06129 

1.57173 

-.16 

-.06930 

1.58866 

-.18 

-.07714 

1. 60560 

-.22 

-.09233 

1.63946 

-.27 

-.11046 

1.68179 

-.32 

-.12077 

1. 72412 

Figure  4.-  Relation  between  s  and  s*. 
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Figure  6.-  Mach  number  as  a  function  of 
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Figure  7.-  Hodograph  of  symmetrical  flow  through  a  nozzle  for 

sQ  =  0.l8(X  =  0.0509117)*  Dashed  lines  indicate  characteristics. 


(a)  Streamlines  and  constant -velocity  lines. 

8.-  Constant-velocity  lines  in  physical  plane  for 
s.  =  0.l8(  A.  =  0.0509117). 
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(b)  Characteristics  (dashed  lines)  and  constant-velocity  lines. 

Figure  8.-  Concluded. 


1.1 
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Figure  9.-  Velocity  distribution  on  streamlines  for  X.  =  0.050912 
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Figure  11.  -  Symmetrical  flow  through  a  nozzle  for 


